R=120; % Сопротивление резистора 120 Ом
C=33e-12; % Емкость конденсатора 33 пФ
Cekv=C/2; % Эквивалентная емкость
tau=2*R*Cekv; % Постоянная времени
T=tau/100; % Ввели шаг
t=0:T:4*tau; % Ввели дискретное время от нуля до 3х тау с шагом Т
lt=length(t);

%Выделение памяти
U=nan(1,lt);
U1=nan(1,lt);
k=nan(1,lt);
Uvh1=nan(1,lt);

%Строим Отклик цепи на постоянное вх. возд.
U(1)=0;
Uvh=1;
for k=2:1:lt; % Среднее значение пропущено, значит единица.
    U(k)=(Uvh+tau*U(k-1)/T)/(2+tau/T);
end;
%step
%АЧХ
%f=0:10:10^6;
%K(f)=1/(1+1i*2*pi*f*(2*R*Cekv));
%G(t)=1/sqrt(1+2*pi*(1/t)*(2*R*Cekv));
%Mag=abs(K(t));
%K(k)=U(k)./Uvh;

%Uvh1=heaviside(k)-heaviside(k-2*tau);

for k=2:1:lt; % Среднее значение пропущено, значит единица.
   if k<lt/2
       Uvh1(k)=1;
   else
       Uvh1(k)=0;
   end; 
    U1(k)=(Uvh1(k)+tau*U1(k-1)/T)/(2+tau/T);
end;

figure(1);
subplot(3,1,1);
plot(t*1e9,U);
grid on;
title('Zaryad C');
xlabel('t, ns');
ylabel('U_c, Volt');

subplot(3,1,2);
plot(t*1e9,Uvh1,'r');
grid on;
title('Heviside');
xlabel('t, ns');
ylabel('U, Volt');

subplot(3,1,3);
plot(t*1e9,U1,'g');
grid on;
title('Otklik na Heviside');
xlabel('t, ns');
ylabel('U, Volt');


